SSC01-VI-6 


Multiple  Quantum  Well  Retromodulators  for  Spacecraft- To- Spacecraft 
Laser  Interrogation,  Communication,  and  Navigation 


N.  Glenn  Creamer,  Ph.D. 

Naval  Research  Laboratory,  Code  8230,  4555  Overlook  Avenue,  Washington,  DC  20375 

gcreamer@  space.nrl.navy.mil 

G.  Charmaine  Gilbreath,  Ph.D. 

Naval  Research  Laboratory,  Code  7215 

gilbreath@rsd.nrl.navy.mil 

Timothy  J.  Meehan 

Naval  Research  Laboratory,  Code  8121 

meehan@ncst.nrl.navy.mil 


Mena  Ferraro  Stell 
Research  Support  Instruments 
Lanham,  MD  20706 
ferraro@ncst.nrl.navy.mil 


Michael  J.  Vilcheck 
Naval  Research  Laboratory,  Code  8123 
vilcheck@ncst.nrl.navv.mil 


William  S.  Rabinovich,  Ph.D. 
Naval  Research  Laboratory,  Code  5654 
rabinovich@nrl.navy.mil 


Abstract. 


This  paper  describes  a  novel  concept  for  laser-based  interrogation,  communication,  and  navigation  between  multiple 
spacecraft  platforms  using  a  gimbaled  laser  source  on  a  pursuer  spacecraft  and  a  target  board  populated  with 
retromodulators  (modulating  retroreflectors)  integrated  on  a  host  spacecraft.  The  combined  laser  source  and 
retroreflectors  can  provide  centimeter-level  relative  positioning  between  each  vehicle,  as  well  as  spacecraft-to- 
spacecraft  laser  communication  via  semiconductor-based  Multiple  Quantum  Well  retromodulators.  Additionally, 
strategies  are  developed  for  utilizing  the  target  board  retromodulator  array  to  provide  relative  attitude  between  each 
vehicle.  In  this  scenario,  each  reflector  has  its  own  unique  modulating  code  sequence,  allowing  the  returned  signals 
to  be  discriminated  and  processed  by  the  pursuer  spacecraft  to  determine  the  relative  orientation.  Based  on 
additional  attitude  sensing  capability,  three  classes  of  host  spacecraft  are  considered:  fully-cooperative,  partially- 
cooperative,  and  non-cooperative.  Numerical  simulations  using  a  five-sensor  target  board  demonstrate  the  potential 
of  the  concept,  and  preliminary  test  results  demonstrate  reflector  discrimination  capability. 


Introduction 

Due  to  the  many  obvious  benefits  of  autonomous 
spacecraft-to-spacecraft  interrogation,  communication, 
and  relative  navigation  for  civilian,  commercial,  and 
military  space  missions,  there  has  been  a  significant 


amount  of  research  and  development  on  associated 
relative  sensor  systems  over  the  past  few  years.1 9 
These  systems  typically  implement  radio-frequency 
communication  links  to  utilize  GPS  sensing  for  long- 
range  relative  positioning  and  combinations  of  visual 
and  laser  ranging  for  short-range  proximity  operations. 
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The  relative  GPS  system  places  the  estimated  target 
location  inside  the  error  sphere  of  the  visual  sensor  at 
about  a  100-meter  range  for  easy  sensor  handoff.  The 
visual/laser  system  is  then  utilized  for  all  subsequent 
short-range  operations  such  as  formation  flying, 
inspection,  and  docking.  Of  course,  this  approach 
requires  the  host  (or  target)  spacecraft  to  possess  an  RF 
antenna  and  transmitter  along  with  a  GPS  sensor  to 
transmit  its  navigational  data  to  the  pursuer  spacecraft. 
Additionally,  this  concept  is  not  suitable  for 
geosynchronous  or  deep  space  missions  due  to  lack  of 
GPS  signal,  unless  a  GPS-type  system  is  emulated 
between  each  vehicle.2  Although  the  visual/laser 
systems  provide  sufficiently  accurate  relative  state 
estimation  for  close-in  operations  (i.e.  0  to  30  meters), 
their  mid-  to  long-range  capability  (i.e.  >50  meters)  is 
significantly  cruder,  with  relative  position  and  attitude 
errors  reaching  10  meters  and  10  degrees  at  ranges 
exceeding  80-100  meters. 

In  this  paper  we  develop  a  novel  concept  utilizing  solid- 
state  multiple  quantum  well  (MQW)  retromodulators  to 
provide  spacecraft-to-spacecraft  laser  communication 
and  navigation  (relative  position  and  orientation).  The 
device  enables  compact,  low  power,  and  low  mass 
optical  data  transfer.  Data  transfer  can  be  on  the  order 
of  megabits  per  second,  depending  on  link 
characteristics,  with  relative  navigation  on  the  order  of 
centimeters  in  position  and  a  few  degrees  in  orientation. 
Links  over  ranges  of  kilometers  down  to  tens  of  meters 
are  possible.  For  proximity  operations  of  about  30 
meters  or  less  (docking  missions,  for  example)  this 
concept  can  work  harmoniously  with  a  vision-based 
system  to  provide  relative  communication  and 
navigation  with  minimal  sacrifice  in  host  power  and 
weight. 

MOW  Retromodulators 
Device  Description 

Modulating  retroreflector  (MRR)  systems  using  MQW 
technology  provide  a  low  power,  low  weight,  multi¬ 
functional  solution  to  the  need  to  reduce  parasitic 
payload  requirements  from  conventional 
communications  technologies.  A  modulating 
retroreflector  is  a  solid-state  device  that  allows  optical 
communication  and  ranging  between  two  platforms. 
MQW  shutters  are  particularly  suited  to  these 
applications  because  the  technology  enables  fast  data 
rates,  requires  very  low  drive  powers,  is  lightweight, 
robust,  and  is  not  polarization-sensitive.11,12  (See  the 
Naval  Research  Laboratory  MRR  web  page  at 
http://mrr.nrl.navy.mil.) 
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Implementation  of  such  a  device  requires  that  only  one 
of  the  platforms  have  an  onboard  laser,  telescope,  and 
tracker.  Thus,  the  device  is  well  suited  to  problems  in 
which  one  platform  has  a  large  payload  capacity  and 
can  serve  as  the  interrogator  and  the  other  platform 
does  not.  The  interrogator  illuminates  the  platform 
carrying  the  modulating  retroreflector  with  a  laser 
beam.  The  laser  beam  is  automatically  reflected  back 
with  no  need  for  pointing  or  tracking.  The  reflected 
return  is  modulated  in  an  On  Off  Keying  (OOK)  mode. 

Bi-directional  communications  can  occur  if  a  lower 
data  rate  is  imposed  on  the  interrogation  beam.  The 
modulated  retroreflected  signal  is  then  received  in  a 
burst  communications  mode.  Detectors  on  the 
retroreflector  platform  can  receive  the  transmitted 
signal,  which  may  inform  the  smaller  platform  of 
interrogator  ID,  or  location  details,  etc.  These  detectors 
can  also  receive  photonic  information  which  can  be 
used  in  acquisition  and  tracking.  A  representative 
concept  is  illustrated  in  Figure  1 . 

The  modulator  must  have  several  characteristics  to 
make  a  link  possible.  The  shutter  must  have  a  high 
switching  speed,  low  power  consumption,  large  area, 
wide  field-of-view,  and  high  optical  quality.  In 
addition,  it  must  work  at  wavelengths  where  good  laser 
sources  are  available,  be  radiation-tolerant  (for  space 
applications)  and  rugged.  Semiconductor  multiple 
quantum  well  modulators  are  one  of  the  few 
technologies  that  meet  all  these  requirements.13,14 

These  devices  are  based  upon  the  same  materials 
technology  as  laser  diodes.  They  consist  of  several 
hundred  very  thin  ( ~  10  nm)  layers  of  semiconductor 
material,  such  as  GaAs,  deposited  on  a  large  (7.6  cm 
diameter)  semiconductor  wafer.  Electrically,  they  take 
the  form  of  a  P-I-N  diode.  Optically,  the  thin  layers 
induce  a  sharp  absorption  feature  at  a  wavelength  that  is 
determined  by  the  constituent  materials  and  the 
structure  that  is  grown.  When  the  device  has  a  moderate 
(-15V)  voltage  placed  across  it  in  reverse  bias,  the 
absorption  feature  changes,  both  shifting  to  longer 
wavelengths  and  dropping  in  magnitude.  Thus,  the 
transmission  of  the  device  near  this  absorption  feature 
changes  dramatically  and  can  serve  as  an  on-off  shutter. 
This  switching  capability  is  shown  in  Figure  2  for  an 
InGaAs-based  MQW  modulator  that  was  designed  and 
grown  for  use  in  an  optical  link.  The  device  is  grown 
on  an  n-type  GaAs  wafer  and  is  capped  by  a  p-type 
contact  layer,  thus  forming  a  P-I-N  diode.  It  is  a 
transmissive  modulator  designed  to  work  at  a 
wavelength  of  980  nm,  compatible  with  many  good 
laser  diode  sources. 
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Unlike  liquid  crystal  modulators,  MQW  modulators 
have  very  high  switching  speeds.  Small  devices 
(diameters  of  microns)  have  been  operated  at  speeds  in 
the  tens  of  GHz.  In  practice,  the  speed  is  limited 
primarily  by  the  RC  time  of  the  device.  Thus,  the  large 
area  devices  (on  the  order  of  a  centimeter)  used  for 
retromodulator-based  communications  typically  have 
speeds  between  1  and  10  Mbps.  Higher  speeds  are 
possible,  however,  depending  on  range  and  the 
sophistication  of  the  fabrication  process.  In  practice, 
data  rates  like  these  are  appropriate  for  many  of  the 
sensors  carried  on  the  small  platforms  for  which  these 
devices  are  intended. 

Data  Link  Considerations 


parallel,  the  power  consumption  may  be  comparable  to 
a  monolithic  device,  but  the  modulation  rate  of  the 
smaller  device  will  be  exploited  while  enabling  the 
larger  aperture.  When  the  drive  voltage  waveform  is 
optimized,  the  electrical  power  consumption  of  a  MQW 
modulating  retroreflector  scales  as 

Dmod4V2B2Rs  (2) 

where  Dmod  is  the  diameter  of  the  modulator,  V  is  the 
voltage  applied  to  the  modulator  (fixed  by  the  required 
optical  contrast  ratio),  B  is  the  maximum  data  rate  of 
the  device,  and  Rs  is  the  sheet  resistance  of  the  device. 
An  example  of  a  waveform  recorded  using  a  9- 
segmented  device  illustrating  10  megabits  per  second  is 
shown  in  Figure  3. 


In  practice,  except  for  close  ranges,  the  link  rather  than 
the  modulator  limits  performance  of  a  modulating  retro 
system.  For  a  conventional  corner-cube  modulating 
retroreflector,  MQW  technology  should  allow  data  rates 
in  the  tens  of  megabits  per  second,  depending  on  the 
range  and  the  interrogator  system. 

For  a  diffraction-limited  system  the  optical  power, 
retro-reflected  from  the  small  platform  back  to  the  large 
platform,  scales  as 


P  D  4  ■  D  2  •  T 

laser  retro  rec  atm 


0div  -PZ 


(1) 


Time-Of-Flight  Range  Capability 

Range  measurements  to  each  MRR  within  an  array  are 
determined  from  time-of-flight  (TOF)  techniques, 
similar  to  the  approach  used  by  laser  range  finders  for 
surveying.  The  unique  aspect  of  our  concept  is  the 
ability  to  selectively  range  each  MRR  with  one  all- 
illuminating  laser  beam.  This  is  accomplished  through 
the  modulation  logic.  The  TOF  is  determined  by 
measuring  the  round-trip  time  it  takes  a  pulse  to  travel 
from  the  laser  diode  to  the  reflector  and  back.  From  the 
speed  of  light,  c,  the  range  becomes 

p  =  ^-(TOF*c)  (3) 


where  Plaser  is  the  power  of  the  laser  transmitter  on  the 
large  platform,  Dretro  is  the  diameter  of  the  modulating 
retroreflector  on  the  small  platform,  Drt.L.  is  the  diameter 
of  the  receive  telescope  on  the  large  platform,  Tatm  is 
the  loss  due  to  transmission  through  the  atmosphere, 
0div  is  its  divergence  of  the  transmit  beam,  and  p  is  the 
range  between  the  two  platforms. 

The  strongest  dependencies  are  on  the  range  and  the 
retroreflector  diameter,  both  of  which  scale  as  fourth 
powers.  Retroreflector  links  fall  off  more  strongly  with 
range  than  conventional  links  because  of  their  bi¬ 
directional  nature.  The  strong  dependence  on 
retroreflector  diameter  occurs  because  increasing  the 
size  of  the  retroreflector  both  increases  the  optical 
power  intercepted  and  decreases  the  divergence  of  the 
returned  optical  beam.  The  link  is  very  clearly  a 
compromise  between  a  large  retroreflector  aperture  to 
maximize  the  returned  optical  power  and  a  small 
modulator  to  maintain  data  rate  while  keeping  the 
consumed  electrical  power  low.  This  trade  is  mitigated 
to  some  extent  by  using  segmented  devices  discussed  in 
the  references.  When  all  the  segments  are  driven  in 


Clearly,  range  errors  are  linearly  related  to  random 
noise  and  bias  errors  in  the  TOF  measurements. 
Through  multi-shot  averaging,  random  noise  errors  can 
be  attenuated  by  the  square-root  of  the  number  of 
samples.  The  bias  errors  are  due  to  un-calibrated  stray 
returns  from  MRRs  that  are  set  in  the  low  or  “off’  state. 

Reflector  Code  Sequence  Discrimination 

Reflector  discrimination  is  achieved  through  a  set  of 
matching  filters  tuned  to  each  MRR’s  unique  modulator 
code  sequence,  as  depicted  in  Figure  4.  The  return 
photons  are  collected  on  the  photodetector  and 
converted  to  digital  signals  via  an  analog-to-digital 
converter.  The  digital  signals  are  then  processed 
through  each  of  the  uniquely  tuned  matching  filters  to 
isolate  the  return  signals  from  each  MRR.  These 
signals  are  then  used  to  determine  the  ranges  to  each 
MRR  and  decode  any  return  data  streams. 

Relative  Pose  Estimation  using  MRR  Arrays 
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The  reflected  signals  from  a  distributed  array  of  MRRs 
can  provide  both  relative  position  and  attitude 
information.  Although  a  variety  of  array  geometries 
can  be  utilized,  we  choose  to  study  the  diamond-and- 
one  configuration  depicted  in  the  photograph  of  Figure 
5,  where  four  MRRs  form  the  corners  of  a  square 
diamond  with  side  length  L,  and  a  fifth  MRR  is 
included  as  a  centering  source  whose  signal  intensity  is 
maximized  during  initial  acquisition  and  track.  As 
described  in  the  previous  section,  each  reflected  signal 
is  modulated  at  a  unique  code  sequence  to  allow 
discrimination  and  signal  isolation  between  each  MRR 
in  the  array. 

We  define  relative  pose  P  of  a  host  spacecraft  with 
respect  to  a  pursuer  spacecraft  as  the  combined  state 
vector  consisting  of  relative  position  states  and  relative 
attitude  states.  In  the  following  we  derive  pose 
estimation  logic  for  three  classes  of  host  spacecraft: 

•  A  Fully-Cooperative  Host  that  provides  complete 
three-axis  attitude  knowledge  to  the  pursuer. 
Examples  may  include  one  member  of  a 
formation-flying  constellation,  a  manned  space 
station,  or  an  orbiting  fluid  storage  platform. 

•  A  Partially -Cooperative  Host  that  provides  a 
measured  body-fixed  unit  vector  of  some  reference 
entity  (i.e.  the  local  geomagnetic  field,  the  Sun 
vector,  or  the  Earth  vector)  to  the  pursuer. 
Examples  may  include  a  minimally-instrumented 
orbiting  commodity  bank  or  one  end  of  a  space 
tether. 

•  A  Non-Cooperative  Host  that  provides  no  attitude 
information  to  the  pursuer.  An  example  may  be  a 
non-instrumented,  passively-stabilized  orbiting 
commodity  bank. 

In  each  case  we  assume  the  pursuer  spacecraft 
possesses  complete  attitude  and  orbit  determination  and 
control  capability. 

Fully-Cooperative  Host 

For  this  class  of  spacecraft  the  host  has  complete 
attitude  determination  capability  relative  to  a  known 
reference  frame  and  incorporates  this  information  on 
the  modulated  optical  signal  returned  to  the  pursuer.  In 
this  case  the  MRR  array  of  Figure  5  is  redundant  and 
we  actually  need  only  one  modulated  reflector  to 
provide  the  necessary  information.  Upon  receiving  the 
modulated  signal,  the  pursuer  spacecraft  calculates  the 
range  to  the  MRR  from  equation  (3)  and  the  relative 
position  vector  from  the  laser  gimbal  angles  and 
calculated  range  using  the  relation 


Rp  =  pi p 


(4) 


where  p  is  the  calculated  relative  range,  RP  is  the 
relative  position  vector  expressed  in  the  pursuer  frame, 

ip 

and  l  is  the  associated  relative  unit  vector  expressed 
in  the  pursuer  frame.  This  unit  vector  is  determined 
from  the  laser  gimbal  azimuth  and  elevation  angles. 
The  relative  attitude  is  obtained  from  the  direction 
cosine  matrix  (DCM)  transformation 


c 


H/P 


=  c  CT 

'-'H/r'-'P/R 


(5) 


where  CH/P  is  the  DCM  of  the  host  frame  relative  to 
the  pursuer  frame,  CH/R  is  the  DCM  of  the  host  frame 

relative  to  the  reference  frame,  Cp  /R  is  the  DCM  of 
the  pursuer  frame  relative  to  the  reference  frame,  and  T 
represents  the  matrix  transpose  function.  Utilizing  this 
complete  relative  pose  information  and  its  own  control 
system,  the  pursuer  spacecraft  can  position  and  orient 
itself  relative  to  the  host  as  required  by  mission 
operations. 

Partially-Cooperative  Host 

For  this  class  of  spacecraft  the  host  possesses  one 
sensor  capable  of  measuring  a  body-fixed  unit  entity 

/v  H 

vector  C  in  the  host  frame  of.  for  example,  the  local 
geomagnetic  field  (from  a  magnetometer),  the  local  Sun 
direction  (from  a  Sun  sensor),  or  the  local  Earth 
direction  (from  an  Earth  sensor).  Depending  on  the 
processing  capability  of  the  host  spacecraft,  either  the 
processed  entity  vector  or  the  raw  sensor  output  can  be 
incorporated  on  the  center  MRR  modulated  return 
signal.  Knowing  the  position  of  the  pursuer  spacecraft 
relative  to  the  reference  frame  from  its  orbit 
determination  capability  and  the  position  of  the  host 
spacecraft  relative  to  the  pursuer  spacecraft  from 
equation  (4),  we  can  determine  the  unit  entity  vector 
measurement  expressed  in  the  pursuer  frame  at  the  host 
location  from  the  relation 

e  —  Cp/Re  (6) 

zv  R 

where  t  is  obtained  from  the  position  of  the  host 
spacecraft  in  the  reference  frame  and  the  known 
position/orientation  of  the  reference  vector  being 
measured  (i.e.  the  geomagnetic  field),  and  Cp  /R  is 
known  from  the  pursuer  attitude  determination  system. 
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Knowledge  of  the  measured  entity  vector  C  in  both  the 
pursuer  and  host  frames  is  insufficient  to  determine  the 
complete  relative  attitude  of  the  host  (two  non-colinear 
vectors  must  be  known  in  each  frame  to  uniquely 
determine  the  relative  attitude).  Therefore,  we  utilize 
the  range  information  from  the  remaining  MRRs  in  the 
array  to  provide  an  additional  vector  measurement. 
Referring  to  Figure  5,  we  choose  the  two  host-fixed 

vectors  rAC  and  rBD  to  maximize  baseline  length. 
Using  the  geometry  of  Figure  6  and  the  law  of  cosines, 
it  can  be  shown  that  the  ranges  to  the  MRRs  defining 
the  host  unit  basis  vectors  satisfy  the  relations 


COS  0[ 


COS  02 


(7a) 


(7b) 


where  the  ranges  are  determined  from  the  modulated 
return  signals  of  each  MRR,  02  is  the  angle  between 

l  and  the  host  h  unit  basis  vector,  and  0  2  is  the 

angle  between  l  and  the  host  h2  unit  basis  vector.  It 
is  clear  that  the  effect  of  errors  in  the  measured  ranges 
(i.e.  random  noise)  will  be  attenuated  by  the  inverse  of 
the  baseline  lengths;  hence,  accuracy  is  proportional  to 
the  length  of  the  baselines  defined  in  the  MRR  array. 
After  solving  for  the  direction  cosines,  the  unit  vector 

l  can  be  expressed  in  the  host  frame  components  as 

lH  =  cos0,  h,  +  cos07  h2 

I - 2 - —  ^  (8) 

-■Jl-cos  0,-cos~02h3 


We  now  have  two  unit  vector  measurements 
represented  in  both  the  pursuer  and  host  frames 

yP  f  II  zvP  aH 

(l  ,  l  ,  e  ,  e  ),  and  can  uniquely  determine  the 
relative  quaternion  (and  DCM,  if  desired)  using  the 
Quaternion  Estimator  (QUEST)  method.13  For  this 
specific  application  the  QUEST  algorithm  provides  an 
estimate  of  the  relative  quaternion  by  optimizing  the 
weighted  loss  function 


J  = 


-C  /p 


+  w2  |eH 


C  e1 


(9) 


6 


where  W  x  and  W  2  are  the  measurement  weighting 
parameters  associated  with  the  laser  vector  and  the 
entity  vector,  respectively.  The  weights  are  typically 
chosen  to  be  the  inverse  of  the  measurement  variances. 
It  can  be  shown13  that  the  optimal  quaternion  qopt 
satisfies  the  eigenvalue  problem 

Kqopt  =^maxqopt  (10a) 


K  = 


S-oI3 

ZT 


z 

o 


(10b) 


^  max  =  VWf  +  W2  +  2W[W2U12  (10c) 


U12  =(/H  -e  )(/  •  e  ) 


f 


+  VH  xeH 


r  xe1 


S  =  Wj{H(/p)T  +  /P(/H)T  } 

+  w2|H(eP)T+eP(eH)T] 


(10d) 


(10e) 


Z  =  w/Hx/p  +  w2eHxep  (iof) 
c  =  w,/H -lp  +w2eH  -ep  (lOg) 


where  1 3  represents  the  3x3  identity  matrix.  Hence, 

upon  solving  for  the  maximum  eigenvalue  from 
equation  (10c),  the  optimal  quaternion  can  be 
determined  from  solution  of  equation  (10a),  subject  to 
the  normalization 

qTq  =  i  (ii) 

Solutions  to  the  eigenvalue  problem  for  the  general  case 
of  n  measurement  vectors  is  also  described  in  Reference 
15. 

Similar  to  the  previous  class  of  host  spacecraft,  the 
pursuer  spacecraft  can  utilize  this  complete  relative 
pose  information  and  its  own  control  capability  to 
position  and  orient  itself  relative  to  the  host  as  required 
by  mission  operations. 

Non-Cooperative  Host 

For  this  class  of  spacecraft  the  host  has  no  attitude 
determination  capability  and,  hence,  provides  no 
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additional  information  to  the  pursuer  other  than  the 
reflected,  frequency-discriminated  signal  returns  from 
each  MRR.  Therefore,  in  addition  to  the  relative 
position  vector  derived  from  the  center  MRR,  the  only 
other  discernible  information  from  the  MRR  array  is  the 
angle  of  the  array  plane  relative  to  the  directional  unit 

vector  l  .  This  angle  can  be  determined  from  equations 
(7a)  and  (7b),  yielding 


(t>  =  y-COS  cos2  0,  -  COS2  02) 


(12) 


where  4*  is  the  angle  between  the  h,  -h2  plane  and 

i 


This  limited  amount  of  relative  orientation  may  very 
well  be  sufficient  for  some  formation-flying  or 
inspection  missions  as  well  as  for  mid-  to  long-range 
phases  of  rendezvous  and  capture  (prior  to  terminal- 
phase  vision  sensor  handoff).  However,  under  certain 
circumstances  additional  information  about  the 
orientation  of  the  host  can  be  postulated.  For  example, 
if  it  is  known  that  the  host  is  passively  gravity-gradient 
stabilized  then,  to  within  a  few  degrees,  it  can  be 
assumed  that  the  Earth  nadir  direction  is  along  the 
minor  principal  axis,  providing  an  approximate  unit 
entity  vector  6  expressed  in  the  host  and  pursuer 
frames.  As  another  practical  example,  if  it  is  known 
that  the  host  is  passively  magnetic  stabilized  then,  to 
within  a  few  degrees,  it  can  be  assumed  that  the  host 
spacecraft  magnet  is  aligned  with  the  local  geomagnetic 
field,  once  again  providing  an  approximate  unit  entity 
vector  in  the  host  and  pursuer  frames.  Given  these 
postulated  entity  vectors,  the  complete  relative  attitude 
can  be  estimated  using  the  approach  described  in  the 
previous  section  for  partially-cooperative  hosts. 

System  Error  Sources 

We  categorize  the  potential  system  error  sources  into 
three  classes:  sensor  errors,  modeling  errors,  and 
geometric  errors.  Sensor  errors  include  combined 
random  noise  and  biases  in  the  time-of-flight 
measurements,  entity  vector  measurements,  and  laser 
gimbal  angle  measurements.  Modeling  errors  include 
host  baseline  knowledge  errors  and  slowly-varying 
biases  due  to  assumptions  on  the  entity  vector  motion 
of  non-cooperative  hosts.  Geometric  errors  include 

effects  due  to  orientation  of  the  laser  unit  vector  l 
relative  to  the  entity  vector  measurement  t  (errors 
increase  as  the  two  vectors  approach  co-linearity)  and 
range-sensitive  angular  errors  from  the  gimbaled  laser 


tracking  system.  These  latter  errors  are  attributed  to  the 
fact  that  the  laser  can  only  be  centered  on  a  return 
signal  to  within  approximately  one-fourth  of  the  total 
beam  divergence  required  to  encompass  all  the  MRRs, 
defined  as  the  baseline  length  divided  by  the  range.  For 
example,  given  a  baseline  of  0.5  meters  at  a  range  of  30 
meters,  the  required  beam  divergence  is  17  mrad  and 
the  laser  pointing  error  is  0.24  degrees.  We’ve 
attempted  to  quantify  these  various  errors  in  Table  1, 
with  the  understanding  that  the  actual  numbers  are 
mission  and  sensor  dependent. 

Kinematic  Filter  Formulation 


To  improve  the  accuracy  of  the  instantaneous  relative 
attitude  estimation  algorithms  given  in  the  previous 
section,  we  develop  a  kinematics-based  Kalman  filter  to 
attenuate  sensor  noise  effects  and  propagate  the  state 
estimate  between  sensor  updates.  We  assume  that  the 
master  spacecraft  has  rate  gyros  to  sense  kinematic 
rotational  changes  due  to  thruster  and/or  reaction  wheel 
torques. 

A  filter  architecture  similar  to  the  Uncoupled  Lefferts, 
Markley,  Shuster  (ULMS)  algorithm16  is  developed 
here.  We  begin  by  defining  the  rate  vector  of  the  host 
relative  to  the  pursuer,  expressed  in  the  host  frame,  as 

(0  =  CO*  -  Co“  =  C0„  -CH/p(0p  (13) 


where  (oj|  is  the  inertial  rate  of  the  host  expressed  in 
the  host  frame,  (Dp  is  the  inertial  rate  of  the  pursuer 

P 

expressed  in  the  host  frame,  and  (Dp  is  the  inertial  rate 
of  the  pursuer  expressed  in  the  pursuer  frame  (as 
measured  by  the  rate  gyros).  An  expression  for  the  rate 
error  vector  relative  to  the  true  rate  vector  (0  can  now 
be  written  as 

(0e  =  co-co 

=  (C0h  -CH/p(Op)-(c0“  -CH(pC0p)  (14) 

=  Aco“  -ACH/p(0p 

where  and  CH  ip  are  the  current  best  estimates  of 

the  host  inertial  body  rates  and  the  relative  DCM, 
respectively. 

The  four-state  quaternion  describing  the  relative 
attitude  of  the  host  to  the  pursuer  can  be  written  as 
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q  = 

q4f 

(15a) 

=0 

(21b) 

q  = 

=  ■&!  q2  q3)T 

(15b) 

#=  FP  +PFT  +  Q 

(21c) 

where  C|  is  the  4x1  quaternion,  q  is  the  3x1  vector 

portion,  and  q4  is  the  scalar  portion.  This  quaternion 
satisfies  the  kinematic  relation 


Q.(  co)  = 


=  y^(CO)  q 


0 

co3 

-co2 

CO, 

—  C03 

0 

CO, 

CO¬ 

co2 

—  (0j 

0 

CO: 

-co. 

-co2 

-C02 

0 

(16a) 


(16b) 


We  can  define  an  error  quaternion  qe  relative  to  the 
true  quaternion  q  from 


q  =  q@qe 


(17) 


0  *I3 

0  0 


(2  Id) 


where  P  and  Q  represent  the  6x6  state  covariance 
and  user-defined  process  noise  covariance  matrices, 
respectively. 

An  approximate  expression  for  the  measured  quaternion 
can  be  written  as 

(q)meas  ®  q  ®  (M)  (22) 


where  (q)  meas  is  obtained  from  the  techniques  in  the 

previous  section,  and  n  represents  a  3x1  noise  vector 
modeled  as  a  zero-mean  Gaussian  process.  Combining 
equations  (14)  and  (19)  leads  to  the  desired 
measurement  model 


where  q  is  the  current  best  estimate  of  q  and  ®  is 
the  quaternion  multiplication  operator.  Following  the 
procedure  in  Reference  16,  it  can  be  shown  that,  to  first 
order  in  the  error  terms,  the  kinematic  relation  for  the 
error  quaternion  becomes 

»  t(A(0h  -ACH/pCO^)-fflxqe  (is) 


(qe)meaS=q"  ®(q)meas 

«  qe  ®  (n,l) 

(23a) 

<<t  l..,,  =  Te  +  n  =  n  | 

%1  +  n 

AcoJJ 

(23b) 

H  =  Jj  0] 

(23c) 

which  can  be  furthered  simplified  by  assuming  small 
rates  for  each  spacecraft,  resulting  in 

4  =  t(A(0h)  (19) 

By  assuming  the  host  body  rates  are  constant,  the 
approximate  error  quaternion  kinematics  can  be 
augmented  with  host  rate  error  kinematic  equations  of 
the  form 

A(8£  «  0  (20) 

The  resulting  Kalman  filter  propagation  equations  for 
the  estimated  quaternion  and  host  rates  take  on  the  form 

yH(ffl)  q  (21a) 


where  (  q"e  )  meas  is  the  vector  portion  of  the  measured 

quaternion  error  obtained  from  equation  (23a)  and  H 
is  the  3x6  measurement  distribution  matrix.  The 
Kalman  filter  measurement  update  equations  can  now 


be  written  as 

K  =PH(HPHT +R)  1  (24a) 

P  =P-KHP  (24b) 

q+  =  q  ®  (qe+ ,  1)  (24d) 


C0»+  =  C&H  +AC0h+  (24e) 
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where  K  is  the  6x3  Kalman  gain  matrix  and  R  is  the 
3x3  measurement  noise  covariance  matrix.  Upon  close 
inspection,  it  is  clear  from  equations  (21)  through  (24) 

that  proper  ordering  of  the  six  states  (<Je,Ad)H) 

uncouples  the  filter  about  each  axis  as  long  as  Q  and 

R  are  diagonal,  resulting  in  three  simple  two-state 
filters  (one  per  axis).  Additionally,  further 
simplification  can  be  achieved  by  implementing  a 
steady-state  or  hybrid  version  of  the  two-state  filters.16 

Numerical  Simulation 

Relative  navigation  of  a  two-spacecraft  leading/trailing 
formation  is  simulated  in  this  section.  The  lead 
spacecraft  serves  as  the  host,  containing  a  target  board 
similar  to  the  one  in  Figure  5  with  one-meter  baselines. 
The  trailing  spacecraft  serves  as  the  pursuer,  whose 
objective  is  to  determine  the  relative  position  and 
attitude  of  the  host  using  a  gimbaled  laser  tracker.  The 
spacecraft  are  placed  in  the  same  orbit  plane  with  an 
altitude  of  600  km  and  an  inclination  of  60  degrees. 
The  host  spacecraft  is  stationed  100  meters  ahead  of  the 
pursuer  spacecraft,  and  both  vehicles  are  aligned  with 
the  local-vertical,  local-horizontal  orbiting  reference 
frame. 

For  this  example  we  assume  the  host  spacecraft 
possesses  either  an  Earth  sensor  to  measure  the  local 
nadir  direction  or  a  magnetometer  to  measure  the  local 
geomagnetic  field.  The  Earth  sensor  has  a  1 -sigma 
noise  error  of  0.1  degrees  and  the  magnetometer  has  a 
1-sigma  noise  error  of  10  nanoTesla.  We  also  assume 
the  1 -sigma  ranging  error  to  each  MRR  on  the  host 
target  board  is  1  cm. 

Three-axis  root-sum-squared  relative  navigation  errors 
are  provided  in  Figures  7  through  9  for  a  Kalman  filter 
solution  using  a  10  Hz  propagation  rate  and  a  1  Hz 
measurement  update  rate.  Unfiltered  position  errors 
with  a  1 -sigma  error  of  about  1  cm,  depicted  in  Figure 

7,  are  consistent  with  the  MRR  ranging  noise.  Both 
unfiltered  and  filtered  attitude  errors,  depicted  in  Figure 

8,  demonstrate  the  noise  attenuation  capability  of  the 
Kalman  filter.  In  the  upper  plot  the  Earth  sensor  is 
utilized,  resulting  in  a  1 -sigma  error  of  about  0.7 
degrees  for  the  unfiltered  solution  (the  gray  plot)  and 
about  0.3  degrees  for  the  filtered  solution  (the  black 
plot).  In  the  lower  plot  the  magnetometer  sensor  is 
utilized,  resulting  in  a  1 -sigma  error  of  about  1  degree 
for  the  unfiltered  solution  and  about  0.4  degrees  for  the 
filtered  solution.  The  twice-per-orbit  error  peaks  and 
valleys  in  the  lower  plot  are  due  to  the  time-varying 
angle  between  the  laser  vector  and  the  geomagnetic 


field  vector,  ranging  from  90  degrees  to  as  low  as  37 
degrees.  In  contrast,  for  the  upper  plot  the  angle 
between  the  laser  vector  and  the  Earth  nadir  direction  is 
consistently  90  degrees.  The  1 -sigma  host  rate  errors  of 
about  0.025  deg/sec  using  the  Earth  sensor  and  0.033 
deg/sec  using  the  magnetometer,  depicted  in  Figure  9, 
demonstrate  the  potential  of  the  Kalman  filter  for 
estimating  host  inertial  rates.  It  should  be  stated  here 
that  only  a  subset  of  the  expected  error  terms  listed  in 
Table  1  were  included  in  this  simulation.  Terms  related 
to  sensor  biases,  modeling  errors,  and  laser  gimbal 
angle  errors  were  not  included. 

Preliminary  Test  Results 

To  demonstrate  reflector  discrimination  capability, 
preliminary  indoor  laboratory  tests  were  performed 
using  a  single  MRR  at  a  range  of  39  meters.17  To 
emulate  two  reflectors,  two  separate  bit  streams  (one 
the  mirror  image  of  the  other)  were  modulated  from  the 
single  MRR,  alternating  at  10  second  intervals.  As 
shown  in  Figure  10,  the  pursuer  transmitter/receiver 
system  consisted  of  a  two-axis  gimbaled  lOOmW  laser 
diode  operating  at  a  976  nm  wavelength,  a  silicon 
photodiode,  signal  amplifier,  analog-to-digital 
converter,  and  two  digital  matched  filters  correlated  to 
each  bit  stream  sequence.  The  host  MRR,  shown  in 
Figure  11,  consisted  of  a  0.5  mm  MQW  retromodulator 
with  a  30-degree  field-of-view,  and  a  three -point 
holder.  The  total  MRR  mass,  including  holder,  was  10 
grams  and  drew  about  75  mW  of  power.  A  15  volt 
swing  was  required  to  achieve  sufficient  on/off  states. 
Reference  17  provides  further  details  on  the  test  setup. 

Signal  acquisition  was  achieved  by  performing  a  series 
of  decreasing  rectangular  searches  based  on  pre-defined 
threshold  signal  levels.  As  the  laser  beam  crosses  the 
MRR  the  return  signal  level  increased  and  a  threshold 
level  was  reached,  initiating  a  smaller  search  pattern. 
This  process  was  continued  until  the  maximum  signal 
level  was  obtained. 

When  a  signal  was  received  by  the  photodetector  the 
receive  algorithm  processed  the  bit  stream  using  the 
matched  filters  to  determine  which  bit  stream  was  being 
sent  by  the  modulator  (in  effect,  discriminating  between 
the  two  streams).  The  outputs  from  the  two  matched 
filters  are  provided  in  Figure  12,  where  it  is  observed 
that  the  on/off  modulation  levels  from  each  bit  stream 
were  properly  matched  and,  hence,  discriminated. 

Future  laboratory  tests  will  utilize  the  target  board  in 
Figure  5,  requiring  five  matched  filters  for  each  MRR. 
In  those  tests  the  MRR  return  signals  will  be  sequenced 
at  a  rate  of  about  0.2  seconds,  requiring  about  1  second 
to  process  and  discriminate  between  each  reflector 
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signal.  That  data  will  then  be  used  to  calculate  the 
relative  pose  of  the  host  platform,  as  described  in  this 
paper. 

Conclusions 

In  this  paper  we  described  a  novel  concept  for 
combined  spacecraft-to-spacecraft  laser  interrogation, 
communication,  and  navigation  utilizing  a  gimbaled 
laser  source  on  the  pursuer  spacecraft  and  a  target  board 
of  modulating  retroreflectors  on  the  host  spacecraft. 
Using  time-of-flight  range  measurements  to  each 
retromodulator  and  the  known  laser  gimbal  angles,  we 
discussed  approaches  for  determining  relative  position 
and  attitude  of  fully-cooperative,  partially-cooperative, 
and  non-cooperative  hosts.  Simulations  demonstrated 
the  potential  of  the  concept  for  providing  centimeter- 
level  relative  positioning  and  degree-level  relative 
attitude  for  vehicle  separation  distances  of  tens  of 
meters  to  kilometers.  Preliminary  experimental  results 
demonstrated  the  capability  of  the  pursuer  detection 
electronics  and  software  for  discriminating  between 
multiple  modulated  retroreflector  signals. 
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Figure  1:  Concept  for  MQW  Retromodulators  used  as  ID  Tags 
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Figure  2:  Transmission  Performance  for  an  InGaAs-Based  MQW  Modulator 
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Figure  3:  Waveform  for  a  9-Segmented  Modulator 
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Figure  4:  Platform-to-Platform  Communication  and  MRR  Discrimination  Strategy 
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Figure  5:  MRR  Target  Board  Layout 
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Figure  6:  Planar  Ranging  Geometry 
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Table  1:  Potential  System  Errors 


Error  Source 

Expected  Value  (3-sigma) 

Sensor  Errors 

Time-of-Flight 

1  -  2  cm 

Entity  Vector  Measurement 

0.1  -  1  deg 

Laser  Gimbal  Angles 

5-10  arcsec 

Modeling  Errors 

Baseline  Knowledge 

1-2  mm 

Entity  Vector  Orientation  (Non-Cooperative  Hosts) 

1-10  deg 

Geometric  Errors 

Co-Linearity  of  Laser  and  Entity  Vectors 

Inversely  Proportional  to  Included  Angle 

Range-Sensitive  Laser  Gimbal  Angles 

Baseline  Length  /  (4*Range) 

4 
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Figure  7:  Simulated  Three-Axis  Root-Sum-Squared  Position  Errors 
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Figure  10:  Optical  Transmitter/Receiver  System  Figure  11:  Mounted  Modulating  Retroreflector 
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Figure  12:  Experimental  Results  Demonstrating  Dual  Bit  Stream  Tracking 
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